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Abstract — In a previous work, the addition of basic screened Coulombic electrostatic 
forces to an existing commercial discrete element modeling (DEM) software was reported. 
Triboelectric experiments were performed to charge glass spheres rolling on inclined planes 
of various materials. Charge generation constants and the Q/m ratios for the test materials 
were calculated from the experimental data and compared to the simulation output of the 
DEM software. In this paper, we will discuss new values of the charge generation constants 
calculated from improved experimental proceedures and data. Also, planned work to include 
dielectrophoretic, Van der Waals forces, and advanced mechanical forces into the software 
will be discussed. 
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I. Introduction 

Theoretical and experimental work was presented on the incorporation of screened Cou- 
lombic electrostatic force into a discrete element modeling (DEM) software used to 
model particle flow [1]. Experimental values of a triboelectric charging constant, a, was 
used to calculate simulated values of Q/m for comparison to experimental values of Q/m. 
The experimental procedure was changed to insure better discharging of the glass 
spheres used prior to the start of each experiment. The resulting values of a give values 
of Q/m in the DEM simulation that compare better to experimental ones. 

Work is proceeding on incorporating advanced electrostatic forces into the DEM soft- 
ware such as the dielectrophoretic and Van der Waals forces. This work is part of a lar- 
ger project with Glenn Research Center (GRC) and DEM Solutions, Inc. to expand the 
capability of the currrent EDEM™ software for planetary as well as terrestrial applica- 
tions. The Kennedy Space Center (KSC) is working on methods to mitigate dust on the 
Moon and Mars [2]. One of these methods uses applied AC electric fields to move the 
dipolar dust simulants off of surfaces. Accurate computer modeling of this electrostatic 
dust screen is a goal of this project. 

II. Revised Experiments 


A. Procedure Modifications 

The experiment consists of an inclined plane apparatus to triboelectrify 2.0 mm diameter 
borosilicate glass spheres that roll down the plane [1]. In the previous experiments, the 
means of neutralizing any initial charge on the spheres was inefficient due to the fact that 
spheres interior to the pile inside the sphere cup were shielded from the deionizer by the 
outer spheres. This resulted in simulation values of Q/m that did not match experimental 
values well [1]. This appartus is shown in Fig. 1. 



Fig. 1. Triboelectric apparatus with aluminum inclined plane, sphere cup, and dam. Total 
charge of the glass spheres, Q , is measured in the Faraday cup shown. 
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To reduce the initial triboelectric charge on the glass spheres, the experimental proce- 
dures were modified to allow better deionization of the spheres prior to each experiment. 
25 spheres were used instead of 500 as before. This number of spheres was easier to 
deionize in the sphere cup prior to the experiment while still providing sufficient spheres 
to observe the flow pattern on the inclined plane. 

The inclined plane was set at 30°. The inclined plane materials used were Aluminum, 
polyvinyl chloride (PVC), low density polyethylene (LDPE), Nylon 6/6, and polytetra- 
fluoroethylene (PTFE). To remove any initial electrostatic charge, the spheres are dis- 
charged by use of an air ionizer. The spheres are then allowed to roll down the inclined 
plane into the Faraday cup by removing the sphere dam. The time of triboelectric contact 
with the inclined plane was measured by a stop watch. The sphere cup, dam, and in- 
clined plane are all of the same material to insure that the triboelectric charging is caused 
by contact only against the test material. 

B. Data 

Ten experiments were performed for each inclined plane material and the results aver- 
aged. The values of Q/m and a from these experiments are summarized in Table 1. 

Table 1 : Values of Q/m and a for revised 30° Experiments (ranges are ± one 

STANDARD DEVIATION ) 


Inclined Plane Material 

Q/m (nCg 1 ) 

« (s' 1 ) 

Aluminum 

0.524 ±0.363 

0.023 ±0.015 

PVC 

5.874 ± 1.610 

0.237 ± 0.042 

LDPE 

0.524 ±0.187 

0.018 ±0.005 

Nylon 6/6 

1.459 ±0.411 

0.057 ±0.012 

PTFE 

2.469 ±0.561 

0.100 ± 0.015 


C. Data Analysis 

Comparison of previous values of a [1] with the values in Table 1 are shown in Fig. 2. 
The much larger a values in Fig. 2 lead to better agreement between the EDEM simula- 
tion and experimental values of Q/m . This data is shown in Fig. 3. These revised ex- 
periments highlight the sensitivity of a to the initial conditions of the surfaces in qustion 
and the physical environment. More on this will be discussed in Sections III A and III B. 
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Fig. 2. Comparison of new and old values of the charge generation constant a. 



Fig. 3. Comparison of EDEM simulation and experimental values of Q/m using the new 
charge generation constant. Q/m values are given in nC g 1 . 

III. Theoretical Aspects 
A. Charge Dissipation Constant 

The data analysis in section II highlights one of the inherent difficulties in using experi- 
mental values of the charge generation constant in that experimentally derived values of 
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the constant can be influenced by initial charge on the surfaces of the contacting bodies. 
This can be especially true when the model equation used assumes zero initial charge on 
the contacting surfaces[l][3] when there is some initial charge. 

Potential parameters of a include, but are not necessarily limited to, temperature, pres- 
sure, % relative humidity, perpendicular force, parallel velocity, surface area, work func- 
tion, and number of occupiable surface sites. The first three are environmental factors, 
the next two are kinematic factors, and the last three are material specific factors. The 
environmental factors can play significant roles in determining the value of a especially 
for insulators [4][5]. More on the initial charge will be discussed in section B. 

B. Revised Triboelectric Charge Generation Equation 

Many tribocharging situations of interest will have initial charge on the particles. To 
account for initial charge on the spheres, Greason’s triboelectric charge generation equa- 
tion [3], 


^- = a{q s -q)-fiq, ( 1 ) 

at 

is re-integrated but this time with limits on the charge integral of q = q 0 at t = 0 in- 
stead of q = 0 at t - 0. Here a and ft are the charge generation and dissipation constants 
respectively and q s is the saturation charge (here taken to be atmospheric breakdown). 
Performing this integration gives 


?(<)=« 



( 2 ) 


The second term of Eq. (3) is the Greason equation disscussed in [1] while the first term 
deals with initial charge on the spheres. Assuming negligible charge dissipation during 
the time frame of the inclined plane experiment (J3 » 0), we get 

q(t)=q 0 e- a ' +q s (l-e al ). ( 3 ) 


For many particles, we multiply Eq. (3) through by the total number of particles, N, and 
rearrange to get 


Nq{t) =Q{t) = Nq s - {i Nq s - Q 0 ]e~ at . (4) 

Here, Q 0 is the total initial charge on all the particles. To use Eq. (4) one would have to 
measure the initial charge on the particles, Q 0 , before the start of any tribocharging ex- 
periment or process. For our inclined plane experiment, we would need to measure the 
total charge on the glass spheres in the sphere cup prior to removal of the dam. A way to 
do this is to replace the sphere cup with a second Faraday cup as the sphere reservoir. 
This will also allow Q/m measurements to be taken both before and after the triboelectric 
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charging. Eq. (4) can be solved for the triboelectric charging constant a. 

t l Nq s -Q 0 _ 


( 5 ) 


Here the parameters N, Q(t), Q 0i and t are measured quantities while the saturation 
charge, q s , for the inclined plane experiment is the electrical breakdown limit of of air 
under standard conditions (2.66 x 1(T 5 C * m" 2 ) [1]. 

C. Electrostatic Force Between a Particle and a Plane 

Another factor to be considered is the electrostatic attraction between the spheres and 
the inclined plane (or other planar surfaces). Initial spheres from the sphere cup will 
encounter an uncharged plane (the plane was discharged prior to each experiment). 
However, the initial rolling spheres will tribocharge the plane and subsequent spheres 
will encounter an oppositely charged plane. The resulting attractive electrostatic force 
will have to be taken into account as a retarding force on the spheres like air resistance 
and friction. 

The electrostatic attraction between a particle and a plane can be physically modeled 
by a charged particle of radius b above a charged circular surface area of radius a. To 
make a point charge approximation, let b« a and the surface charge density on the cir- 
cular area, a, be uniform and constant. This is shown in Fig. 4. 



Figure 4. Charged particle of radius b near a charged circular plane of radius a. 
In cylindrical coordinates, the unit vector R is given by [6] [7] 
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R - —rr + z 0 z 


with 


R _ -rr + z 0 z 
*| yR +z o 


( 6 ) 


Using Eq. (6) and the geometry from Fig. 4, the differential force element can then be 
written as [6] 


dF = 


4 7T£ 0 R 


-rr + z 0 z 


i 


r 2 +zl j 


prdrd<f > . 


(7) 


Before Eq. (7) is integrated, it should be noted that all the radial components will can- 
cel from symmetry. Integrating Eq. (7) gives 


F = 


qa_ 
4 ne n 


ft 




( 


\ 


( r2+z o) 


,3/2 


drd</> z = 


2 s n 


1 

1 


+ z l J 


z. (8) 


Letting Zo » b for close proximity between the sphere and other surfaces we get 

f , A 


F = 


qa_ 

2^n 


1- 


V 


4a*+b 2 


z . 


(9) 


Eq. (9) can be used to approximate the electrostatic force between a charged particle 
and a charged surface. 


IV. DEM Simulation 

DEM Solutions, Inc. has performed additional particle-particle and particle-surface tribo- 
electric simulations modeling a physical situation similar to the experiment described in 
Section II. In this simulation [12], four angles (0°, 15 °, 30 °, 45 °) are modeled simulta- 
neously. This geometery is shown in Fig. 5. Box structures are used as the particle fac- 
tories. The generic particles (20 for each factory) modeled are 0.5 mm in diameter with a 
mean density of 1200 kg * m' 3 and a shear modulus of 7.7 x 10 6 Pa. The charge genera- 
tion constant, a, was set at -0.5 s' 1 to insure that the particles would charge negatively 
and the sloped trough would charge positively. 
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Figure 5. EDEM Model of four inclined curved surfaces at 0°, 15°, 30°, and 45°. Box 
structures on the end are the particle factories. 

Charge build-up on the inclined surfaces are calculated in EDEM and are shown graphi- 
cally in Fig. 6. The higher charge build-up on the lower inclines is due to the longer con- 
tact time of the spheres. Additional normal force due to gravity in the lower inclines is 
negligible with such small particles but could be a factor for larger ones. 




Figure 6. Charge build-up on the inclined surfaces at 0°, 15°, 30°, and 45°. 
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The zero degree slope generated the highest charge even though relatively little particle 
motion occurred [12]. Average charge on the spheres for each angle is shown in Fig. 7. 



Figure 7. Average charge on the spheres for each angle [12]. 

The charge on the spheres at zero degrees in Fig. 7 is unphysical and is a consequence 
of contact time between the spheres and the inclined surfaces, not triboelectric charge 
exchange. This highlights a limitation of Greason’s charge generation equation [3] in 
that it assumes rolling or sliding triboelectric contact between the two surfaces. The 
EDEM code will have to be revised to take non-motion contact into account. 

V. Future Work 


A. Additional Experiments 

The experiments described in Section II need to be revised further to insure the removal 
or measurement of any initial charge on the spheres. This will include changing the ex- 
perimental proceedure to discharge the spheres as they are being loaded into the sphere 
cup by a deionizer and to load them slowly. Also, as mentioned in Section II B, the 
sphere cup can be replaced with a Faraday cup to measure the initial Q/m of the spheres 
prior to release and the resulting data compared to the previous experiments. 

Additional materials will be used for the experiment to increase the data set. Addi- 
tional sphere materials are: Copper, Carbon Steel, Acrylic, Polytetrafluroetylene (PTFE), 
High Density Polyethylene (HDPE), Carbon, and Nylon 6/6. These spheres range in 
diameter from 1 .6 mm to 2.34 mm. 

The number of spheres will be reduced from previous experiments [1] from 500 to 
250. This number will still allow flow patterns to be observed on the plane and facilitate 
the discharge of the spheres prior to the start of the experiment. 
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B. Addition of the Dielectrophoretic and Van der Waals Forces 

KSC and Glenn Research Center (GRC) are engaged with DEM Solutions, Inc. to add 
advanced electrostatic and mechanical forces to the DEM software produced by DEM 
Solutions, Inc. KSC is assisting DEM Solutions, Inc. to add the dielectrophoretic force 
to the software. The dielectrophoretic force is the force exerted on the induced dipole 
moment of an uncharged dielectric particle by the presence of a non-uniform electric 
field [8]. This is the force that moves the dipolar Lunar and Martian regolith simulants 
off of KSC’s dust screen [2], The basic equation of the dielectrophoretic force is 


F DEP — 2.7T£\R 


s 2 -e x 




+ 2 


\VE 2 


( 10 ) 


where £j is the electric permittivity of the medium, e 2 is the electric permittivity of the 
particle, R is the particle radius, and E is the non-uniform electric field [8]. Eq. (10) is 
valid for most non-uniform electric fields but if the non-uniformity of the field is very 
large (an approximate definition for “large” is where the non-uniform electric field is 
stronger than other fields or forces present), then a multipolar treatment is necessary 
[8][9]. The multi-polar force acting on a sphere with the electric field in the z-direction 
[8] is given by 


F_ = 


2 xe x KR 2n+ ' 

8 


dz 


d n ~'E, 


dz 


n - 1 


( 11 ) 


where K = — — 

s 2 -i- 2e x 


(also called the Clausius-Mossotti factor). 


When n = 1, Eq. (10) is recovered. The question to be answered is how intense of a 
field gradient is required for higher multipolar terms of Eq. (11) to become significant? 
Reference [9] indicates that higher multipoles of the dielectrophoretic force are not neg- 
ligible within a few particle diameters unless the external field is weak. 

Van der Waals forces is the term generally ascribed to intermolecular forces other than 
covalent or electrostatic [1 0][ 1 1 ]. In more macroscopic systems, sub-micron particles 
can adhere to each other or to surfaces via Van der Waals forces alone. GRC will pro- 
vide the necessary theoretical support to DEM Solutions, to incorporate these forces into 
EDEM and provide experimental verification of the resulting simulation data. 
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